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The shell model solve the nuclear many-body problem in a restricted model space and takes into 
account the restricted nature of the space by using effective interactions and operators. In this paper 
two different methods for generating the effective interactions are considered. One is based on a 
partial solution of the Schrodinger equation (Bloch-Horowitz or the Feshbach projection formalism) 
and other on linear algebra (Lee-Suzuki). The two methods are derived in a parallel manner so that 
the difference and similarities become apparent. The connections with the renormalization group 
are also pointed out. 
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I. INTRODUCTION 

The shell model has been in use for many years. 
There are two main variants. The first is due to 
Bloch and Horowitz This method is closely re- 
lated to Feshbach projection operator formalism^] 
frequently used to the justify the optical model. 
The renormalization group treatment of the two- 
body scattering due to Birse et al 3] is also closely 
connected to this method. The second method is 
due to Lee and Suzuki p|. It is used in the no-core 
shell modeljs] and to derive the nucleon-nucleon 
effective interaction V\ ow & [E ■ It is the purpose 
the present paper to put the differences and sim- 
ilarities between the two methods into sharp con- 
trast by a parallel derivation of the two methods. 

II. THE OPERATOR u> AND ITS USES 

We start with a projection operator P = P 2 and 
its compliment Q = 1—P. Following Lee-SuzukiQ 
we define an operator: 

QuP\^) = Q\*P) (1) 

While this is a common part of the Lee-Suzuki 
method it is rarely used in the Bloch-Horowitz 
method. Here we use it in both methods and as 
will be seen the main difference between the two 
methods is precisely in how u> is determined. First 
we note that only the off-diagonal matrix elements 



of lu are needed. In general u> is not well defined — 
a projection operator does not have an inverse. We 
need independent information on for example 
that it is a solution of the Schrodinger equation 
or a member of a well define set. However once u> 
is known, we can easily construct effective inter- 
actions and operators. For example, starting with 
the Schrodinger equation H\tp) = E\yb) we have: 

H{P + QujP)\i/j)=E\ub) (2) 

PH(P + QuP)\tp) =EP\i/)) (3) 
(P + PjQ)H(P + QuP)\i/;) 

= E(P + PjQ)\ij) (4) 

= E(P + P^QuP)\ip) (5) 

where eq.n nas been used to eliminate Q\ip). The 
form, eq. |21 is the simplest and is the one uni- 
versally used in the Bloch-Horowitz method. It 
has the advantage that Li occurs linearly. The last 
has the advantage of being symmetric. The states 
P\tf>) are only orthogonal with with respect to the 
weight function. Orthogonal functions can be ob- 
tained by defining 

\<t>) = y/P + P^QuP\^) . (6) 

The effective interaction is now given as (see 
ref. 0): 

1 = =(P + P^Q)H 
y/P + Puj^QluP 

x(P + Qu;P) I g |0)=£|0).(7) 

\JP + Puj^QojP 
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The equations, [3] through [7| are formally equiv- 
alent. Which one to use is matter of taste and 
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convenience. In particular using \ip) has the ad- 
vantage that as we increase the model space we 
just add additional contributions leaving the terms 
already present unchanged whereas 4> has to be 
renormalized at each stage. On the other hand 
the <p have the advantage that they are orthogonal 
(at least for some choices or lj) and an hermitian 
operator is being used. 

By construction, eq.|3|gives an effective interac- 
tion in the P space that reproduces the P space 
components of the wave function from the full- 
space calculation. If P commutes with the kinetic 
energy (or more generally the unperturbed Hamil- 
tonian) then the half off shell t-matrix is also re- 
produced since: 

T\iPp W ) = PV\i>) = V(P + QuP)\i>) (8) 

where \ippw) is a plane wave. Acting on the left 
side of the last equation with P restricts the t- 
matrix to its low momentum components below 
the cut off. So we have PT\^ PW ) = PV(P + 
QluP)\iP) and can readily identify the low energy 
effective potential that reproduces the half off-shell 
t-matrix as: 



V eff = PV(P + QujP). 



(9) 



This will be true for either the Bloch-Horowitz or 
Lee- Suzuki choice of w. As we will see shortly, the 
Bloch-Horowitz method reproduces the fully off- 
shell t- or k-matrix leading to the same results as 
Birse et aZ's|3j application of the renormalization 
group. On the other hand, in the derivation of the 
effective interaction, Vi ow the authors only re- 
quire the half-off shell t-matrix to be reproduced. 
That approach is equivalent 6] to the Lee-Suzuki 
method. 

It is equally easy to generate the effective oper- 
ator for a given initial operator O: 

<^|0|^> = (il>\(P + PJQ)Q{P + QuP)\i/>) (10) 
1 



y/P + PuJ^QluP 

(P + QloP) 



(p + pjQ)e 



y/P + Puj^QluP 



10} (11) 



As these examples illustrate u> is the key ingredi- 
ent. Once it is known, the effective quantities can 
be calculated easily and the formal connections 
become clear. The great advantage of the present 
method using u> is that many of the derivations 
become very simple. 



III. THE BLOCH-HOROWITZ METHOD 

The first method for solving for u> is act on eq. |2 
with Q, use eq. ^and solve for to. This gives: 



?wP|V> 



1 



E-ib — QHQ 



-QHP\i>) (12) 



There is one aspect of this equation, that thought 
of as an operator equation, is quite odd (although 
so common that one hardly notices it anymore). 
Namely the right hand side depends explicitly on 
the energy of the state — the energy has been 
given a ip subscript to emphasize this point. The 
equation is not a general operator relation but is 
only valid for a limited set of matrix elements — 
those with the given eigen function of H on the 
right hand side. This is, presumably, the source of 
the starting energy problem: What energy do we 
use once we take this equation out of this narrow 
context? In spite of all that, if is an eigen 
function of H then ea. 1 121 must be satisfied. Even 
the alternate form for u) derived in the next section 
must obey this equation. 

Eq-Elis au that is needed to derive the effective 
interaction from eq. |3| Inserting eq. El m to eq. [3] 
we obtain: 



[ PHP + PHQ 



E,p - QHQ - ie 



QHPj P\ij}) 
= EPtPIU) 



This equation can be immediately recognized as 
the Bloch-Horowitz shell-model equation or the 
Feshbach optical model equation. A bit of a mir- 
acle has occurred here since eq. |3] is asymmetric 
with to appearing only on the right hand side while 
eq. 1131 is symmetric. Nothing is gained by going 
to eq. 01 or [3 The effective interaction is not her- 
mitian due to the ie. Since to is energy depen- 
dent both the effective interaction and the effective 
operators are energy dependent. For off-diagonal 
matrix elements different energies and hence w's 
are needed on the left and right hand sides of 
eq. El The energy dependence of the effective 
interaction is also a major contributor to making 
the wave functions non-orthogonal. The effective 
interaction and operators can be used for any state 
with a non-zero overlap with the model space. 
This is in contrast to the Lee-Suzuki method dis- 
cussed in the next section. 

The effective interaction, eq. El can be de- 
rived equally well starting with the Lippmann- 
Schwinger equation T = V + VGqT. Defining 



3 



T P = PTP and T Q 
gives the equation: 



QTP and eliminating Tq 



\PHP + PHQ 
[PHP + PHQ 



1 



- QHQ - le 
1 

En - QHQ - ie 

xG T P 



QHP) 
QHP 



Thus the Bloch-Horowitz effective interaction re- 
produces the low momentum components of the 
fully off-shell t-matrix (or cquivalently the k- 
matrix). For two-body scattering, it reproduces 
the results of Birse et al 3] . 



IV. THE LEE-SUZUKI METHOD 



The second method |4j, |9j] of obtaining w is a bit 
more subtle. We divide the Hilbert space into two 
parts \k) and \k'). For the sake of definiteness take 
the states to be eigenf unctions of the full Hamilto- 
nian. Consider an orthonormal set of states, \olp), 
that span the P space. If the states P\k) form 
a complete basis (not necessarily orthonormal) in 
the P space then we can write any one of the \ap) 
states as: \ap) — ^ fc P\k}a^ F . We now act on 
this equation with {oiq\Qu)P where (q;q| is an ar- 
bitrary state in the Q space. Using eq. Q]we have: 

(a Q \QuP\a P ) = ^2(a Q \QujP\k)a^ 

k 

= Ewx" ( 14 ) 
k 

Since the \otp) are an orthonormal set the expan- 
sion coefficients a? p can be determined by invert- 
ing the equation: 



{a P \ap) 



E 



(a' P \k)a?*> . (15) 



This is the defining equation of the Lee-Suzuki 
method. Once this equation is solved we have ui 
and can calculate the effective interactions and op- 
erators. However the inverse of (a' P \k) must exist. 
If not you have to choose a different set of |fc)'s or 
modify the projection operator P. The existence 
of the inverse implies that there are the same num- 
ber of states in the P space and the set |fc). 

The oddity of the Lee-Suzuki method is that 
ui depends not only on P but also on the states 
chosen to be in the set For example, even if 



you are only interested in the lowest state, the ui 
acting on it will change if the set is modified 
to include the 101st state rather than the 100th. 

Note the dramatic differences compared to the 
previous method. Much to the surprise of people 
use to working with the Bloch-Horowitz or Fcsh- 
bach method the effective operators in this method 
are energy independent. Also this method can 
work for no more states than the number of states 
in the P space. This is also in distinct contrast to 
the Bloch-Horowitz method. 

Since the V\ ow k obtained in the renormaliza- 
tion group method is equivalent to the Lee-Suzuki 
method you might expect that a matrix inversion 
like that defined in eq. 1151 would also play an im- 
portant role in the renormalization group deriva- 
tion. Indeed this is the case, see for example eq. 4 
and 5 in ref. . The renormalization group V\ ow k 
is given by eq. |5|with current choice of ui. 



V. UNITARY TRANSFORMATIONS 

It is well known that the Lee-Suzuki method 
corresponds to a unitary transformation 9] and 
eq.0does suggests a unitary transformation with: 



U = 



1 



v/1 + PuAQujP 



(P + PjQ) 



(16) 



However, U is not a unitary transformation since 
UW — P rather than one. What U does is take 
a state in the \k) space and transform it into the 
P space while W does the inverse, namely maps a 
state in the P space into a state in the |fc) space. 
The operator, U, acting on a state in the \k') space 
gives zero. A full unitary operator would have to 
map states in \k') into the Q space. An operator 
that does this is: 



U = 



1 



y/l + Plu'iQujP 
1 



v/1 + QluPlu^Q 

1 - (QujP - Plo^Q) 



(P + Plo^Q) 
(Q - QujP) (17) 
(18) 



v/1 + Plu^QujP + QluPuj^Q 

Note (QojP-P^Q) 2 = -{P^QujP + QujPuj^Q). 
When acting on a state in the \k) space this re- 
duces identically to U since Q — QujP acting on 
state in |fc) gives zero. Thus if one is only in- 
terested in states in the \k) space then the only 
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difference between the two operators is aesthet- 
ics and ink usage. The full form does indicate 
that what the Lee-Suzuki method is doing, cither 
explicitly or implicitly (depending on the formula- 
tion) is constructing a unitary transformation that 
block diagonalizes the the Hamiltonian. 

For a state in the set \k), eq^Jcan be simplified. 
Using eq. ^ and the definition of cj> (eq. EJ) gives 
U\k) = y/l - Pus^QusPP\k) = \<f> k ). Using the 
completeness of the \(f>k) in the P space \<j>^) can 
be expanded as \4>k) — Yli P\l){k\P\l)~ 1 / 2 . In this 
form the \4>k) are manifestly orthonormal. The 
wave functions P\k) span the P space but are not 
orthonormal. The \4>k) have been orthogonalized 
in a way that treats all the states in the \k) space 
symmetrically. 

In the Bloch-Horowitz method it also possible 
to construct a unitary transformation however it 
is much less useful since us depends on the energy 
of the state of under consideration and the unitary 
transformation would be different for each state. 



VI. CONCLUSIONS 

The important quantity in determining effective 
interactions or operators is us. A projection oper- 
ator loses information. It is the job of us to re- 
store the lost information. The parallel derivation 
of the Bloch-Horowitz and Lee-Suzuki methods 
shows that the essential difference between them 
is how the lost information is recovered or equiv- 
alently how us is determined. The crucial step in 
the Bloch-Horowitz method is a partial solution of 
the Schrodinger equation while in the Lee-Suzuki 
method the key step is a matrix inversion. 



While the main equations 10 through |5J) in the 
two methods are similar there are significant dif- 
ferences due to the different methods of restoring 
the lost information. The Bloch-Horowitz method 
deals with one state at a time and has energy de- 
pendent effective interactions and operators. It 
can describe any state with a non-zero overlap 
with the model space. The Lee-Suzuki method 
deals with a set of wave functions, a set with the 
same dimension as the P space. It can describe 
only those states in that set. The effective inter- 
actions and operators are energy independent. 

Another method for using a restricted space is 
the renormalization group. The crucial equation 
in comparing to the renormalization group is eq.|21 
It shows you how to construct a model-space ef- 
fective operator that gives the same model-space 
components of the wave function as the calcula- 
tion in the full space. Either the Bloch-Horowitz 
or Lee-Suzuki methods can be used. The first re- 
produces the fully off-shell t-matrix and is thus 
equivalent to the procedure of Birse et al. The 
second only reproduces the half-off shell t-matrix 
and is equivalent to V\ ow k- 
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